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We review four facility location problems which are moti- 
vated by urban service applications and which can be thought 


of as extensions of the classic @-median problem on networks. 


In problems P1 and P2 it is assumed that travel times on net- 
work links change over time in a probabilistic way. In P2 

it is further assumed that the facilities (servers) are movable 
so that they can be relocated in response to new network 
travel times. Problems P3 and P4 examine the Q-median 
problem for the case when the service capacity of the facili- 
ties is finite and, consequently, some or all of the facilities 
can be unavailable part of the time. In P3 the facilities have 
stationary home locations but in P4 they have movable loca- 
tions and thus can be relocated to compensate for the unavail- 
ability of the busy facilities. We summarize our main results 
to date on these problems. 


Keywords: Facility location, network, Q@-median problem, 
travel times 


1. Introduction 


The location of facilities is a spatial allocation of 
scarce resources. Cities, whose residents require a 
large number of diverse services, pose particularly 
difficult questions of facility location. Some facilities 
are fixed, such as libraries, schools, and out-patient 
clinics, whereas others are mobile service providers 
such as police vehicles and taxi cabs; still others have 
mobile service providers (or servers) at fixed facilities 
or home locations, such as fire engines and pumpers 
at fire stations and ambuiances at hospitals. Some 
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require each service request or customer to travel to 
the facility; others require the facility or a server 
located at a facility to travel to each customer. 

When attempting to apply in cities classical deter- 
ministic location theory methods such as those 
derived from Q-median and/or Q-center (where Q 
is the number of facilities to be located) formulations, 
one often encounters the following problems: 

(1) Travel times on urban streets are probabilistic 
quantities influenced by traffic conditions, weather, 
street repair work, traffic lights, etc. 

(2) Customer arrival times and locations, as well 
as their associated service requirements, are also prob- 
abilistic quantities, creating various levels of spatio- 
temporal congestion in the system servicing these 
customers. 

(3) Facilities do not have unlimited service capa- 
city and thus may be unavailable at the time a nearby 
service request occurs; this, in turn, creates a pattern 
of workload sharing, in which facilities other than 
the closest facility — which may be temporarily 
unavailable — quite regularly service customers. 

(4) Urban service systems are often multi-objective 
systems in which acceptable system performance is 
determined by a set of measures of efficiency, equity, 
and effectiveness. For instaice, for an ambulance sys- 
tem a measure of efficiency is mean city-wide travel 
time for a given number of ambulances; a measure of 
equity is the magnitude of the inaximum difference 
between neighborhood-specific mean travel times; 

a measure of effectiveness is the number of lives 
saved or hospital days reduced by actions taken by 
the ambulance crew. 

Thus, the classical location formulations are 
inapplicable in complex urban settings because they 
assume a Single-objective, deterministic world in which 
the closest facility is always the one utilized. 

In recent years, we and our colleagues have been 
attempting to incorporate some of the complexities 
cited above into new, usually probabilistic formula- 
tions of multi-facility location problems. Our approach 
has been to attempt to bring into the analysis one or 
two such complexities at a time, recognizing that a 
fully comprehensive model and location algorithm 
remain beyend tractability at this time. Yet, it is 
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Table 1 


Comparison of the Q-median problem with the four problems discussed in this paper 


P1: Medians or 
stochastic networks 


Classic 
Q-median problem 
networks 


Demand at nodes Demand at nodes 


P2: Movabie servers 
on stochastic 


Demand at nodes 


P3: Congested 
medians 


P4: Congested movable 
server systems 


Demand at nodes Demand at nodes 


only only only only only 

Unlimited service Unlimited service Unlimited service Limited service Limited service 
capacity capacity capacity capacity # capacity # 
Fixed-location Fixed-location Movable location Fixed-location Movable location 
facilities facilities facilities or servers # facilities or servers servers 4 

Deterministic Probabilistic Probabilistic Deterministic Deterministic 

link lengths link lengths # link lengths 4 link lengths link lengths 

Service by Service by Service by Service by Service by 

closest facility closest facility closest facility closest available closest available mobile 


or by closest mobile 


server 


facility or by server 4 
closest available 


mobile server 7 


4 Assumption different from that of the classic @-median problem. 


our hope that progress along each of several fronts 
will begin to indicate the types of new location results 
that obtain when further realism is brought to loca- 
tional modeling. 

In this paper, we present a review of results that 
have been obtained recently ona set of four prob- 
lems, labelled P1 through P4 (see Table 1), which in 
many ways are similar to the Q-median problem on 
networks [8,9]. 

For comparative purposes we summarize the Q- 
median problem as follows: 

(1) There are Q facilities to be located on a net- 
work G(N, L). 

(2) Requests for service (‘customers’) occur only 
at the nodes NV. 

(3) Each request is serviced by the closest facility, 
which is assumed to be always available to provide 
service, 

(4) Travel time on any fraction 0 of a link is 
assumed to be equal to (total link travel time) - @ 
(0<@ <1). 

(5) The problem is to find the set of Q facility 
locations that minimizes mean travel time associated 
‘with a random service request. 

We have modified the Q-median problem to 
include such considerations as probabilistic service 
demands and travel times, unavailability (or conges- 


tion) of facilities and mobility of either facilities 
or servers located at facilities. Referring to Table 1, 
the problems may be summarized as follows: 


Pl: Medians on stochastic networks, This is essentially 
the Q-median formulation, but with the condition 
that travel time on each link is a discrete random 
variable. Our formulation, which is precisely detailed 
in Section 2, allows link and thus network travel 
times to change only at discrete time epochs and 
assumes that the traveller (either the customer or the 
server) has perfect knowledge of current network 
travel times. In applications, the discrete time epochs 
may refer to time intervals between receipt of reports 
of traffic conditions. Given any particular set of 
values for link travel times, it is assumed that each 
customer is serviced by the closest (in current travel 
times) facility. The problem is to locate the Q-facili- 
ties in order to minimize mean travel time associated 
with each customer, where the average is taken over 
all network (iravel time) states. P1 has been dealt 
with in detail elsewhere [8,17] and thus will be 
treated here only briefly as a special case of P2. 


P2: Movable servers on stochastic networks, This is 
the same as P1 with the additional complication that 
the location of each facility may be changed as net- 
work travel times change. For instance, during rush 
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hour conditions, the city of Boston’s Department of 
Health and Hospitals may change the location of one 
or more of its ‘downtown’ ambulances if severe traffic 
congestion occurs; due to a tunnel, which may become 
a very congested link during rush hours, one option 
that is currently employed is to move one of the 
downtown ambulances to tne exit side of the tur nel 
(outside of downtown) during such saturation condi- 
tions. Here the discrete time instants for possible 
changes in travel times may represent ‘periodic loca- 
tion review times’ associated with assessing the possi- 
bilities of ambulance relocation during rush hours. 
Model P2 however assumes unlimited service capa- 
city, with service by the closest facility, thus limiting 
its applicability to those ambulance services having 
low ambulance utilizations (i-e., low probabilities 

of being unavailable). Model P2’s assumptions are 
somewhat less restrictive when considering certain 
types of customer-to-facility-systems, including 
mobile outpatient clinics, medical testing services, 
lunch vans, book-mobiles, etc. The optimization 
problem with P2 is to find the optimal set of facility 
locations for each network state and the associated 
relocation rules that direct facilities when the net- 
work changes state. The cost function to be mini- 
mized is a weighted sum of the mean travel time to 
customers and of the cost of facility relocation. 


P3: Congested medians. The formulation described 
here changes focus from probabilistic link lengths to 
probabilistic availability of servers (or facilities). It is 
the same as the Q-median problem with the compli- 
cation that each facility (or server) may or may not 
be available to provide immediate service. Service 

is provided by the closest available server or, if all 
servers are busy, the customer is handled by some 
back-up service system. This model is most relevant 
to emergency service systems such as police, fire, 
emergency medical and emergency repair. The opti- 
mization objective is to find the set of @ facility loca- 
tions that minimizes mean travel time associated with 
a random service request, where the average is com- 
puted over all possible combinations of servers avail- 
able and unavailable. A special case of this formula- 
tion is the hypercube model [14], which has been 
utilized by several cities in North America and Europe. 


P4: Congested movable server systems. This formula- 
tion is similar to P3, with the additional complexity 
that servers are movable and a server that is en route 
to a location may be assigned (dispatched) to a cus- 
tomer. This model combines all of the esseritial ingre- 


dients of actual ambulance systems except probabi- 
listic travel times. The movable servers in this instance 
are allowed in order to be able to respond to tempo- 
rary conditions of congestion in certain parts of the 
network; called relocation or move-up in fire services, 
the aim of relocating servers is to maintain in real 
time at least some minimal level of coverage in each 
area of the city (or network). The possible assign- 
ment of moving servers to customers assumes perfect 
knowledge of real-time server locations. To date, 
emphasis in P4 has been on determining a compute- 
rized server relocation policy for a given set of 0 
possible locations. 


Throughout our paper emphasis is on model 
assumptions and results. Thus, no formal proofs will 
be given, but appropriate references are provided for 
the interested reader. [lustrative examples have also 
been included. 


2. P2: Mobile servers on stochastic networks 


Let GW, L) be an undirected network with N the 
set of nodes ({V| =) and L the set of links. Service 
demands are generated exclusively at the nodes of 
G(N, L) with h; being the conditional probability 
that a demand comes from node i (¢ € NV) given that a 
demand was generated (7; can be viewed as the ‘nor- 
malized weight’ of the node /). 

In problem Pz travel times on the network change 
in a probabilistic manner according to the following 
simple model: At constantly spaced intervals (epochs) 
G(N, L) is assumed to undergo changes of state. If 
rand s are two distinct states of the network and if 
tg(i, /) indicates the travel time on link @, j/)@Z when 
the network is in state 8, then t,(, /) #¢,(i, /) for at 
least one link (i, /)€ L. Let Mf be the set of all possible 
states of G, |M| =m. 

Transitions between network states at the epochs 
are governed by an ergodic Markov transition matrix 
P with p,, © P being the probability of a transition 
from a state r to a states 7 EM, s EM). We also 
denote the steady—state probability vector of the 
matrix P as Il (ILP = I, 7? , I, = 1). The Q mobile 
servers (of infinite service capacity) which are to be 
located on the network are operated as follows: When- 
ever there is a demand for service, that customer is 
assigned by a system operator and travels to the 

server closest to it in terms of travel time. A customer 
not using such a system operator is assumed to know 
the current travel times and server locations, always 
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selecting the closest server. Whenever there is a change 
of state of the network, the operator of the service 
has the option of relocating one or more of the ser- 
vers. A relocation of a server is associated with a cost, 
which we shall choose to express in units of travel 
time. The operator’s objective is to minimize the 
long-term expected cost (again expressed in terms 

of units of travel time) of providing the service. The 
long-term cost per epoch will be a weighted sum of 
the total expected travel time of demands to the 
servers per epoch (under all states of the network) 
and of the expected cost of the server relocations 
that take place per epoch. 

We now define some additional quantities needed 
to express this problem in quantitative terms. Let 
K(r) = {k10), k2(7), «... Kg(r)} be a set of Q points, 
where the @ servers are located whe: the network is 
in state r. We shall denote the shortest travel time 
between any two points x and y on G when the net- 
work is in state 7 as d,(x, y); the shortest travel time 
between any point in the set K(r) and a specific 
point x on G when the network is in state s as d,(K(r), 
x); and the shortest travel time between the ath 
point in the set K(r) and the yth point in the set 
K(s) (for a and y = 1, 2, ..., @) when the network is 
in state / as d,(K,,(r), K,,(s)). 

The cost (in units of travel time) of relocating the 
ath server in K(r) to the yth location in K(s) with the 
network in state s is given by f[d,(K,(r), Ky(s))]. We 
also define binary variables W,(K4(r), K(s)) as fol- 
lows: if the server at K,(r) is relocated to the location 
K,(s) when the state of the network changes from r 
to s, then W,(K g(r), K,(s)) = 1; otherwise it is equal 
to 0. ; 

Finally, we define as a strategy, any vector K = 
(K(1), K(2), .... K(m)) with m elements, where each 
element K(r), r © M, provides the set of Q locations 
where the servers will be placed when the network 
is in state 7. A simple strategy is any strategy with 
K(1) = K(2)= + = K(m), ie., a strategy in which 
servers remain stationary under all states of the 
network. 

We shall now state the assumptions under which 
the results for this model have been derived: 

(1) The travel times ¢,@, /) for all ry © M and all 
a@A)EL are finite. 

(2) The time required to travel a fraction 0 (0 <6 
< 1) of any link (@, 7) EL fe: any state r EM is 
equal to @ - ¢,(i, /). 

(3) The current state of the network is known to 
the service operator (or, in the absence of an opera- 


tor, to the customers) at all times. 

(4) Time intervals between changes of state are 
muci: longer than trip times on the network. 

(5) No demands or further state changes occur 
while servers are being relocated after a change of 
network state. 

(6) All servers are available whenever a demand 
occurs. (This assumption is the same as that in the 
classic Q-median problem.) 

(7) The relocaticn cost function, f( - ), is non- 
decreasing concave. 

(8) Demands are always assigned to the closest 
mobile server. 

Assumption | asures connectivity of the network 
under all states. A less restrictive version of Assump- 
tion | (one that allows some of the ¢,(i, /) to be infi- 
nite while still leading to the same results) is given in 
[1,16]. Assumption 2 concerning uniformity of 
travel time on any given link is necessary in the proof 
of Theorem 1. Assumption 3 allows all travel in the 
network to be along the shortest path. Assumption 4 
renders negligible the probability that link iravel 
times will change while a customer is travelling to a 
server. (Were this to happen the server might no 
longer be travelling on a shortest travel time path.) 

Assumptions 5 and 6 are the major simplifying 
assumpiions in P2, Both assumptions would be approx- 
imately true, in practice, if, for instance, the average 
time interval between generation of demands on the 
network was much longer than travel times on the 
network (assuming that demands are generated accord- 
ing to a stationary renewal process independently 
of the state of the servers). As mentioned in the Intro- 
duction, Assumption 6 would be approximately true 
if each movable server in a customer-to-server system 
could simultaneously serve any number of customers 
(e.g., mobile libraries or ‘bloodmobiles’). 

Assumption 7 is necessary for Theorem { to hold. 
It implies ‘economies of scale’ for the cost of travel 
times — a reasonable hypothesis in most practical 
contexts. (Obviously the family of acceptable func- 
tions, f( - ), also includes the linear cost function.) 
Assumption 8, stating that the closest server to a cus- 
tomer always services that customer, is identical to 
that of the standard Q-median problem. 


2.1, The problem 
We can now express our objective function in terms 


of :he quantities that have been defined. For any 
given strategy K = (K(L), K(2), .... K(v)), the quantity 
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mt n 
A=2y Ly 1,hd{K(n), i) (1) 
r=) 1 
gives the long term (‘steady—state’) expected travel 
time per dispatch or server assignment. Similarly the 
quantity 


m om 2 @Q2 
Bec Qi 2s Mypalls Le WK), Ky) 
r=1 1-1 a=] y¥=1 
l#r 
Xfldf{Ka(r), Ky()I | (2) 


represents the long-term expected cost of server relo- 
catioris per transition epoch, taking into account all 
possible changes of state from any possible state. The 
constant c in (2) is the relative weight assigned to the 
expected server relocation costs per epoch, in order to 
express relocation cosis in units of travel time. 

If an average of \ customers arrive per transition 

—3zpoch, our problem is to minimize 


Z=ATB. (3) 
2.2. Main result 


Our main results for problem P2.can now be sum- 
marized, beginning with the following Hakimi like 


Theorem 1. At least one set of optimal locations for 
P2 exists on the nodes of the network. 


The proof of Theorem 1, although long, is quite 
straightforward (see [1]). Basically it uses the proper- 


State | 


ties of concave functions to show that Z cannot 
increase if a server which is located on some link 

(i, 7) = L£ when G is in some state r € M, is moved to 
one of the nodes / or /, no matter where the other 

Q — 1 servers are located. Theorem 1 reduces the total 
number of strategies to be considered to (G)”, since 
it reduces the number of candidate locations to (9) 
for each of the m states of the network. It also leads 
to a straightforward integer programming formulation 
of P2 [1]. The size of this integer programming prob- 
lem grows very quickly as n, Q, and m increase. We 
have little computational experience with problems 
of this type to date. However, due to considerable 
similarity with the formulation of the deterministic 
median problem, the recent research of Erlenkotter 
[5], Garfinkel, Neebe and Rao [7], Revelle and Swain 
[18], Galvao [6], Jarvinen, Rajala and Sinervo [11], 
and Cornuejols, Fisher and Nemhauser [4] can per- 
haps be of help in solving our problem as well. 


2.3. An example 


The following simple example illustrates some of 
the ideas above. The network of Fig. 1 can be in one 
of the two states, 1 and 2. The numbers next to the 
links represent lengths (travel times), whereas those 
next to the nodes are the weights, h;. Obviously, the 
only difference between the two states is in the travel 
time on link (3, 2). Customers appear at a rate of 
A = 1 customer per transition epoch. 

The Markovian transition matrix, P, that describes 
the statistical dependence between the two states is 
shown below. 


Fig. 1. The network un-ler states 1 and 2. 
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The steady state probabilities associated with P are 
Il, = 0.4 and I, = 0.6. 

Let f(y) = by be the relocation cost function, 
an increasing concave function (b is a positive con- 
stant) and let c = 1. Suppose that we wish to locate 
two servers on the network. It is easy to show that 
for 0 <b <0.2357 the optimal strategy is {K(1), 
K(2)} = {G, 4), 1, 4)}, ie., one of the servers must 
be relocated from node 3 to node 1 when the net- 
work makes a transition from state 1 to state 2 (and 
vice versa when the change of states is from 2 to 1), 
while the second server remains stationary at node 4 
irrespective of the state of the network. Note that 
by moving from node 3 to node 1 (and back) as the 
network changes states, the first server retains some 
proximity to node 2. 

On the other hand, when b > 0.2357 the optimal 
strategy is a simple one {K(1), K(2)} = {(1, 4), (1, 4}. 
Note that, in this case, demands at node 2 are serviced 
by the server at node 1 when the state of the network 
is either 1 or 2 but not via the same shortest path. 


2.4, Additional results 


Ongoing research has led to some additional find- 
ings regarding P2 [3]: 

(1) For the case in which a single mobile server is 
to be located on a probabilistic tree, the optimal strat- 
egy is to keep the server stationary at a single node, 
independently of the state of the tree. That node is 
the (single) median of the tree which remains unchang- 
ed under ail states of the tree. 

(2) For the case in which a single mobile server is 
to be Iecated on a probabilistic network with m 
states, we have devised a heuristic algorithm which, 
in effect, amounts to sclving a succession of classical 
1-median problems, i.e. 1-median problems on deter- 
ministic networks. These latter problems can be 
solved efficiently and in a straightforward way (see, 
for instance, [8]). 

(3) Simple upper and lower bounds on the mini- 
mum value of Z can be found for the most general 
case of Q mobile servers with m network states. By 
sc.ving m independent Q-median problems, one for 
each network state, both an upper and a lower bound 


on Z can be obtained. A second upper bound is the 
value of Z associated with the best available simple 
strategy (this is also the solution to P1, a special case 
of P2). 


3. P3: Congested medians 


For P3 and for the remainder of this paper, travel 
times on the links of G(W, L) are assumed determinis- 
tic. (This also permits some simplification in our nota- 
tion.) The concern in P3 is with locating Q stationary 
facilities each having one server which can become 
busy (or unavailable). The objective is to locate the 
Q facilities so as to minimize the expected travel time 
associated with a random service request weighted 
appropriately by the equilibrium (steady-state) state 
probabilities of the system. In contrast to P2, ‘states’ 
here are defined according to the status (busy/avail- 
able) of each of the facilities (and has nothing to do 
with travel times on the links). 


3.1, A one-server example 


To motivate our discussion of P3 and to demon- 
strate some general characteristics of spe ially distri- 
buted queuing systems, we consider first a simple 
one-server example: 

Suppose that an ambulance is to be garaged at a 
hospital located at point x on a roadway (link) con- 
necting towns (nodes) a and b (Fig. 2). A fraction 
hg (hy) of ambulance demand is generated from town 
a (b). (tg + hy = 1.) No demands occur on the link. 
In response to each demand, or service request, the 
ambulance, if available, travels to the appropriate 
node at a speed v, spends a fixed time 7 at the scene, 
then travels back to the hospital (with the patient on 
board), again at a speed v. If the ambulance is not 
immediately available, due to servicing an earlier 
demand, then the current service request is entered 
into a queue of other delayed service requests. 


Location of ambulance and hospital 


0) x as 


Fig. 2. Single link, two-node ambulance example. 
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Queued requests are serviced in a first-come, first- 
served manner. The service-time, S, associated with a 
service request is the sum of the travel time to the 
scene, the on-scene time, and the travel time back to 
the hospital. 

The mean and variance of the service time are 
readily computed, 


E[S] =S =7+ [1 +2x(a, — hy) /v, (4a) 
E[(S — S)?] = 02 = 4x?[1 —(tg — hy)? |v? (4b) 


Now suppose that the total demand rate, summed 
over both towns, is A ambulance requests per unit 
time. The ambulance utilization factor, or fraction of 
time the ambulance is busy, is 


p=aS =dfrt [1 + 2x(hq — hy) /v} , (5) 


In order for an equilibrium (steady-state) solution to 
exist, we must have p < 1. Note that certain locations 
x, for given d, v, hg, and hy, may violate this assump- 
tion, causing the queue to grow without bound; where- 
as other locations may satisfy p <1. 

In this example, we want to locate the hospital 
and ambulance at a point x* € [4 ; 4] which minimi- 
zes the mean total response time (in-queue plus 
travel time) required for the ambulance to reach a 
random service request. if Wg is the mean in-queue 
time, then this mean response time can be written 


Ta = Wq + 5 [1 + 2x(hg — hy)I /v- (6) 


Now to complete our specification cf this loca- 
tional problem, we need only to describe the nature 
of the stochastic process that generates service requests 
from a and b. This, in turn, determines W,. Motivat- 
ed by physical considerations, we assume that each 
town generates ambulance demands independently 
as a homogeneous Poisson process at a rate Mg 
(town a) or Af, (town b) per unit time. The com- 
pleted model is a single-server queue having Poisson 
input and general independent service times, abbre- 
viated M/G/1 in queuing parlance. For such a queue 
it is well known that [13] 


pS +o? 
So 7 
q 2(1 — p) ( ) 


Using (5) and (7) with (6), we can obtain an expres- 
sion for the mean response time Tp as a function of 
given system parameters and the locational decision 
variable x. Our objective is to select x* € [-4,4 
which minimizes (6). 


It is a straightforward exercise in differential cal- 
culus to find the required minimum. The essential 
point is that for certain sets of values of system 
parameters the optimal location is at one of the 
nodes; for other sets of values, however, the optimal 
location is at some interior point x © (-4, 4). For 
instance, for a symmetric system in which hg = hy, = 
4 the mean service time does not depend on x; thus 
p does not depend on x. The only quantity in (6) 
that depends on x is the quadratic in the numerator 
of Wy, 02 = 4x?/v*, which is clearly minimized at 
x = 0. Thus, assuming p < 1, for a symmetric system 
(i.e., hg = h,) of the type we are discussing, the opti- 
mal facility location is always at the link mid-point, 
x = 0. Moreover, this optimal location moves away 
from x = 0 continuously as (4, — h,) moves contin- 
uously from 0. 

To illustrate typical behavior, consider an example 
for which 7 = 1, v= 1, A= 4. If the ambulance were 
located at x = 0, we would have 


S=r+i/l=1+1=2; 02 =0;p=}-2=4. 


(Note that one-half of the mean service time in this 
instance is on-scene time and the other half is travel 
time). Regardless of hg and h,, an ambulance located 
at x = 0 would provide the following service perfor- 
mance: mean travel time to the scene = 4; Wg = 
[4 - 2]/24) = 1; 7p =1+4=2. Thus, any move- 
ment away from x = 0 (for hg #h,) must reduce 
Tr below 3, Such movement would increase o2 
above 0, but may decrease S and p, resulting in a net 
decrease of the objective function Tp. The optimal 
location x* is shown as a function of (A, — hg) in 
Fig. 3. Note that for |k, — 4, | 2 0.3, the optimal 
location is at the busier of the two nodes; for all 
other values of (hy — hg), x” is an interior point. In 
fact, if no queuing ever occurred, the objective func- 
tion to emerge is simply 4[1 + 2x(hg — hy) /v, a 
linear function whose minimum exists at the busier 
of the two nodes a or b (or if hg = hy, at any point 
x & [-4,4]). This result also obtains if a queue is 
tot allowed to form, by assigning customers arriving 
when the server is busy to a back-up system, where 
‘cost’ per response is assumed higher than that of the 
primary system. 

The observations from the M/G/1 analysis have 
motivated cur general treatment of problem P3, 
having multiple-servers. Our resv!ts to date include 
only a class of subproblems for wiich nodal solutions 
are found to be optimal. However, the reader must 
not think that general spatially distributed queuing 
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Fig. 3. M/G/1 example: optimal location of facility as a func- 
tion of hy — hg. 


systems will always exhibit locational optimality 
at nodes. 


3.2. General formulation 


Let Xq be the set of all possible locations of Q 
facilities (Q > 1), on the network G, i.e., 


Xg = {Xg = 1, XQ) X, EG, K=1, aa Qh}. 


Given any location Xg = (x), .... xg) © Xq, let X, 
denote that the facility at x, is not staffed with an 
available server (the facility is busy) and %;, that the 
facility x, does have an available server. Therefore, for 
any Xg € Xo there are 22 combinations (states), any 
one of which the network can assume, according to 
the status of the @ facilities. Let Y x(q) be the set of 
all states for Xg € Xg and let y xg) (or for convenience 
yg) be a generic element of Y xg). 

Let ¢(i, /) be the travel time on link (i, /), @ JEL, 
and let d(g, j) be the (minimum) travel time asso- 
ciated with the closest available server ' to node j, 
when the system is in state yg. 

Our results for P3 have been obtained under the 


1 Such a zero-lookahead strategy is very reasonable, but not 
always optimal in the sense of minimizing time-average mean 


following set of assumptions: 

(1) Demands (service requests) occur according to 
a general renewal process. 

(2) Each request requiies a total service time (=on- 
scene service time plus travel time) whose distribution 
is general and not dependent on the identity of the 
server or the history of the system. This implies that 
on-scene service times are assumed to be much grea- 
ter than travel times so that variations in total service 
times which are due solely to variations in travel 
times among potential servers can be ignored. 

(3) Appropriate ergodicity conditions apply so 
that a unique steady state distribution exists. 

(4) The time required to travel a fraction 0 (0 <6 
< 1) of any link (i, j) EL is equal to @ - t(i, /). 

Assumption 2 is reasonable for many urban ser- 
vice systems which are often characterized by on- 
scene service roughly an order of magnitude greater 
than typical travel times. Assumption 4 is used in 
the proof of Theorem 2 below. 

For any possible set of locations XY g © XQ, let 
P(yqg) be the steady-state probability that the net- 
work is in state Yyx © Y x(qy. Let yQ be the state in 
which all the Q facilities are busy (i.e., vO = (1, X2. 
.. ¥g) in our notation). 

Conditioned on any state yg € Yx¢gy ~ {yh}, 
the expression 


n 
Qu hjd(vo, i) 
q= 


is the expected travel time associated with a random 
service request. If, however, the system is in state 
YO (i.e., all Q servers are busy), when a service 
request occurs, we will assume that the request is han- 
died by a back-up system having a travel time cost 
R per response (R assumed larger than any travel 
time in the primary system). (We have also obtained 
results for systems that allow queuing, under specific 
queuing service disciplines and assuming negative 
exponential service times.) 

The congested median problem is now stated: 


Xx gex Q 


travel time. An optimal policy occasionally requires assign- 
ment of other than the closest available server [12], in or- 
der to leave the system in a state which best anticipates 
future service requests. We do not consider such strategies 
in our formulation of the congested median problem. 
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with 


F(Xg)= Pog) 2a hyd(vo.1) 


a 
yoEY x(Q)- 9} 
+ P’9) -R. 


Obviously, the classic median problem is a special 
case of (8) arising when P(g) = 0, Vyvg #1, «+s 
xg) — the state where all the units are evailable. The 
weights P(y’g) in (8) represent the fraction of time 
that the network is in each of the 22 possible states. 
Therefore, as noted before, we take into account the 
fact that any subset of facilities can become depleted 
of servers. 


3.3, Results 


Our main results for P3 begin once more with a 
Hakimi-like theorem: 


Theorem 2. At least one optimal solution to (8) 
exists on the nodes of the network. 


The proof of Theorem 2 is given in [2]. We believe 
that it can be extended to more general dispatch poli- 
cies, using criteria other than ‘closest-available server’ 
to assign demands to servers. 

Once the solution space has been limited to the 
nodes of G, the following important observation is 
useful in solving the congested median problem for 
specific situations: Whenever 

(a) demands over the nodes of G are generated in 
a Poisson manner at a total rate A and at each nodej 
independently in a Poisson manner with rate A; = AA, 

(b) the iotal service time for each server is nega- 
tive exponential with mean 1/u, and 

(c) each facility can hold exactly one server, then 
(8) can be viewed as a locational objective function 
for the hypercube queuing model [14], (see [15] for 
an overview of the hypercube model). 

In [12] Jarvis developed a heuristic algorithm for 
finding a set of, hopefully, near optimum locations 

in the framework of the hypercube model where loca- 
tions are constrained to nodes and each node can: 
contain not more than one facility. 

For the purposes of P3, it should be noted that, 
because of Theorem 2, both the hypercube model 
and Jarvis’ algorithm do not suffer from a loss of 
generality by considering locations only on nodes. 


We have applied Jarvis’ algorithm, in the manner des- 
cribed by Larson [15], to several moderate-size net- 
work problems (e.g., 7 = 25, Q = 4) using (8) as the 
performance measure to be optimized. Convergence 
of the algorithm has usually been fast and the results 
have been surprising on occasion. It is not uncommon 
to find a solution consisting of a set of locations which 
is among the weakest possible sets in terms of the clas- 
sic (deterministic) median problem. 


3.4, Example 


The following example will illustrate some of our 
previous discussion. Suppose we want to locate three 
facilities on the simple network shown in Fig. 4. The 
numbers next to the nodes are the fractions of 
demands from each node /, / = 1, ..., 5 and the numbers 
next to the links are the travel times. There are (3) 
possible distinct locations: 


{1,2,5}, {1,3,4}, {1,3, 5}, 
{2,3,5}, {2,4, 5}, {3, 4, 5}. 


Adgdi sis. leeahh 
{1,4,5}, {2,3, 4}, 


The optimal location according to the standard 3-me- 
dian problem is {1, 2, 5}, which can be obtained by 
hand computation. Suppose however that service 
requests occur in the network in a Poisson fashion 
with A = 4, and the service time for each one of the 
three units is exponential with identical means u~! = 1. 
Let us assume a zero capacity queue with R = 5 units 
of time — the cost resulting when dispatching the 
reserve unit. We also assume that server preferences 
are determined solely by geographical proximity. 

The Jarvis algorithm with an initial location at the 
absolute 3-median, i.e., {1, 2, 5}, converges after one 
iteration to the optimal solution at location {2, 3, 5}. 
The improvement achieved by moving from the loca- 
tion {1, 2,5} to {2, 3, 5} is 3% in terms of the objec- 
tive function of the congested median problem. It is 
interesting to realize that the location {2, 3, 5} is 
among the weakest possible locations in terms of the 


Fig. 4. A simple 5 node network. 
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standard median problem. This indicates that blind 
application of the classic (deterministic) median prob- 
lem can lead to erroneous results, even for such 
simple networks. 


4. P4: Congested systems with mobile servers 


This problem was the main topic of research of 
one of the authors’ doctoral dissertation (see [1, 
Chapters 3, 4, and 6]). It is also the problem (among 
P1—P4) with which our computational experience 
is most extensive. Due to space limitations our dis- 
cussion of P4 will be limited and mostly descriptive. 

In P4 our concern is not so much with initially 
locating Q mobile servers but in exploring how free 
servers at any given time should be repositioned to 
‘compensate’ for servers which are busy and, there- 
fore, unavaiiable. The objective is to minimize a 
weighted sun of the expected travel time to a ran- 
dom requesi for service and the expected server re- 
positioning costs per unit of time. It is always assumed 
that each request for service is serviced by the - 
closest available server, with che exception of 
requests which find all servers busy. These latter 
requests are presumed lost and a penalty, R, is paid 
whenever such an event occurs in a manner analogous 
to that described earlier for problem P3. A related 
assumption is that an idle server is eligible for being 
assigned to a request for service, even while in the 
process of moving, from one location to another. 

In order to provide a proper focus for indicating 
the type of approach taken as well as the limitations 
of our research, let us consider a case in which Q = 3. 
Assume that somehow it has been determined that 
the best set of locations for the three (indistinguish- 
able) servers, whenever all three are available, is at 
the nodes {i, j, k} and that at a particular moment _ 
the three-server system has just entered state (i, j,k) 
as described by the status of each of the three servers, 
in the notation described under P3. If servers are 
restricted to locations i, j, and K only (these, for 
instance, might be the locations of fire houses) then 
the following three actions are possible when the sys- 
tem enters state (i, 7, kK): 

(i) do nothing; 
(ii) move an idle server from node j to node i; 

(iii) move an idle server from node k to node i. 

If action (i) is taken, the possible successor states 
to (i,j, k) are (i,/, &) if the next event to take place 
is the completion of service to the call with which 


the server from / is currently occupied and (i, j, k) or 
(,7, 4) if the next event is a request for service 

closet to j or to k, respectively. The immediate cost 
associated with ‘do nothing’ is the sum of the expect- 
ed travel times associated with servicing demands 
from / or from k (as the case may be) weighted by 
the probabilities of transitions to (¢ j, k) or to 

(é, j, k) respectively. 

If action (ii) is taken the possible successor states 
are: (i, j, k), if the next event is a completion of tke 
busy unit’s service; (i, 7, k) if the next event is a 
request for service and the server which was moving 
from j to i happened (in mid-trip) to be closer than 
the available server at kK to the request for service, at 
the instant when the request for service occurred; 
(i,7, k) if the next event is a request for service and 
the server at K is closer to the request for service than 
the moving (from / to 7) server at the instant when the 
request for service occurs. An important assumption 
is implicit in the statement above: during the reposi- 
tioning period only one completion of service or arri- 
val of a new service request (but not both) is possible. 
A direct consequence is that if a service request occurs 
during the repositioning period, and if the moving 
server is not dispatched to it, then the next service 
request can occur only after the repositioned server 
reaches its destination. This assumption facilitates 
the analysis that focuses on time instants following 
a transition. The expected immediate cost for action 
(ii) takes into consideration the possibility that the 
moving server will be dispatched to the next service 
request. The situation for action (iii) is entirely symume- 
trical to that for action (ii). 

The last two paragraphs discussed transitions, 
transition probabilities and immediate costs associ- 
ated with state (/,7/, &). A similar analysis can be per- 
formed for all seven other states. Note that for states 
(i 7, k) and (i,j, k) only the ‘do nothing’ action is 
available. , 

We are now in a position to present our assump- 
tions and outline the results that we have obtained 
for P4: 

(1) Demands over the nodes of G are assumed gen- 
erated in a Poisson manner at a total rate A and at 
each node j independently in a Poisson manner with 
rate Aj = Aji. 

(2) The total service time for each server is nega- 
tive exponential with mean 1/y. 

(3) 1/u >> (travel times on the network) so that 
all variations in service times due to the travel times 
can be ignored. As noted under P3 this assumption is 
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reasonable for several important urban services. 

(4) All servers travel at the same constant speed 
on the network G. 

(5) Only one server can be repositioned at any one 
time. This implies that, for @ = 4, when two servers 
are busy and the other two idle, it is not permitted 
to reposition simultaneously the two idle servers to 
the positions of the two busy servers. (This, however, 
does not preclude repositioning first the one server 
and then the other.) 

(6) Perfect information about the position of 
any moving (due to repositioning) or stationary 


servers is available at all times to the server dispatcher. 


Assumption 6 is reasonable in view of the avail- 
ability of increasingly accurate and reliable ‘automat- 
ic vehicle location’ (AVL) systems for urban service 
systems. Unfortunately, the only reason for assump- 
tions 4 and 5 is that we have found the analysis to 
be mathematically intractable without them. 


4.1. Results 


As the discussion of states, state-transition proba- 
bilities, alternative actions and immediate costs may 
have suggested, we have been able to approach and 
solve the version of P4 described above as a Marko- 
vian decision process (MDP). For any given set of Q 
nodai locations {f;, i2, ....ig} housing the Q servers 
(the ‘home locations’), this solution finds the opti- 
mum repositioning policy that minimizes a weighted 
sum of the expected travel time to a random request 
for service and the expected server repositioning 
cost per unit of time. This objective function is for 
MDP’s the ‘average cost per unit of time’, usually 
denoted as g. An optimum policy is a listing of the 
actions to take for each possible state of the system 
in order to minimize g. 

A rather involved computer program (about 900 
APL statements) has been written to accomplish 
this. The program calculates all transition probabili- 
ties and immediate costs for all possible states and 
permissible actions and then applies Howard’s policy 
iteration algorithm [10] to find the optimum policy 
and the associated minimum value of g. We have 
solved many problems of moderate size (nm = 25, Q = 
3) in order to explore the sensitivity of the optimum 
policy and of g to the server utilization ratio p(=)/Qu) 
and the “goodness” of the home locations !. Reduc- 


2 Since we are dealing with a zero-line capacity queue, p = 
A/Qu is not the average function of time servers are busy. 
This latter quantity is less than p due to lost calls. 


tions of up to 20% in the objective function, g, over 
the ‘all-do-nothing’ policy (i.e., no repositioning 
allowed) have been observed in the examples attempt- 
ed. The greatest benefits from repositioning have 
been observed at the middle-range of server utiliza- 
tion rates, which is also the most interesting and 
common range for the applications with which we 
are concerned. The reasons for this are quite obvious: 
at very low server utilization rates, dispatching to 
requests for service is so infrequent that the bene- 
fits resulting from reduced travel times due to server 
repositioning are too small to justify the effort of 
repositioning; by contrast, at very high utilization 
rates, the servers are busy most of the time, leaving 
little room for maneuverability by the operator of 
the service. In our numerical examples we have also 
found cases in which the optimum policy is unex- 
pected and “‘counter-intuitive” (see Section 4.2). 

An interesting by-product of this research has 
been the development of an efficient algorithm for 
keeping track of which available server is closest to 
each node of G as an available server moves (due to 
repositioning) between any two nodes on G. 

We have found that the analysis of P4 with Q = 2 
is relatively simple and that some special simplifying 
results can be derived for this case. In fact, for Q = 2, 
we have been able to develop results for the case in 
which the expected service times are different for 
each of the two-servers and for the case in which, 
when one of the two servers is busy, the other one 
can be repositioned to any node on G and not only 
to the other available home location. 


4.2. Example 


In this section we present an example intended to 
demonstrate that for high utilization factors and/or 
for large travel times along the repositioning paths, 
repositioning often occurs in anticipation that the 
moving unit will be near ‘strong’ nodes while travel- 
ling and will be dispatched to the next service request. 
To show that let us refer to the network in Fig. 5. 

Let us consider the case where the two existing ser- 
vice units are located at nodes 1 and 24 and R = 40. 
It is easy to verify that node 24 is ‘stronger’ than 
node 1 (225, hj d(24, /) = 15.84 < 72, hy d(1,/)= 
19.10). 

The model was applied for several cases when the 
utilization factor p varies from zero to infinity. As 
expected, for low utilizations the optimal policy 
is to move the service unit from node 1 to 24 when 
the network enters state (1, 24) and to do nothing 
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Fig. 5. An urban transportation network. 


when the network enters state (1, 24). For high uti- 
lization, however, the optimal policy is to reposition 
the available server for both states (1, 24) and (1, 24). 
To understand this result, note first that the distance 
between those two nodes is quite large — 30 units of 
time. In addition to that, all the nodes within a dis- 
tance of less than 26 units from node 24 are ‘stronger 
than node 24. In fact, the shortest path form node 24 
to node 1 includes the two strongest node: node 10 
and node 14 (which is also the absolute median). 
Therefore if the decision is to reposition for state 

(i A 24), then unless p is very small, chances are high 
that the repositioned server will be dispatched from 

a ‘stronger’ position than node 24. Only for very 
small utilizations is it likely that the repositioned 
server from node 24 might be dispatched near or at 
node 1, in spite of the large distance between nodes 

1 and 24. 


> 


5. Conclusions 


We have reviewed in this paper a class of problems 
which begin to bridge the gap between deterministic 
location formulations as represented by the median 


problem and an important area of urban applications 
characterized by probabilistic demands and travel 
times, facility unavailabilities (congestion) and server 
mobility. Several potentially useful results have been 
obtained under sometimes restrictive sets of assump- 
tions. By examining the entries in Table 1, represen- 
ting the types of complexities we have added to the 
standard Q-median formulations one can readily 

see that more work is needed in integrating the types 
of methods used in this paper. For instance, a model 
incorporating both probabilistic travel times and prob- 
abilistic availabilities of servers would be useful. Even 
for the complexities analyzed here, more work is 
needed on reducing the restrictiveness of certain 
assumptions and on devising computationally efficient 
algorithms. Finally, future work should also begin to 
include the all-too-important multi-objective nature 
of many urban service systems. 

We are hopeful that the continued analysis of 
these problems will tie more closely together pre- 
viously disparate efforts on network analysis, queuing 
theory, spatial analysis, and Markovian decision pro- 
cesses. 


References 


[1] O. Berman, Dynamic repositioning of mobile servers on 
networks, Technical Report No. 144, M.LT. Operations 
Research Center, Cambridge, MA (1978). 

(2] O. Berman and R.C. Larson, The congested median prob- 
lem, Working Paper OR 076-78, M.I.T. Operations 
Research Center, Cambridge, MA (1978). 

13] O. Berman and A.R. Odoni, Locating mobile servers on 
a network with Markovian properties, Working Papcr 
OR 083-78, M.I.T. Operations Research Center, Cam- 
bridge, MA (1978); revised version. 

{4] G. Cornuejols, M.L. Fisher and G.L. Nemhauser, Loca- 
tion of bank accounts to optimize float: an analytic 
study of exact and approxiinate algorithms, Manage- 
ment Sci. 23 (1977) 789-810. 

{5] D. Ertenkotter, Dual-based procedure for uncapacitated 
facility location, Working Paper 261, Western Manage- 
ment Science Institute, University of California, Los 
Angeies, CA (1976). 

[6] R.D, Galvao, A dual-bounded algorithm for the ?- 
median problem, Paper presented at the International 
Symposium on Locational Decisions at Banff, Alberta, 
April 24-28, 1978. 

[7] R.S. Garfinkel, A.W. Neebe and M.R. Rao, An algorithm 
for the m-median plant location problem, Transporta- 
tion Sci. 8 (1974) 217-236. 

{8] G.Y. Handle: and P.B. Mirchandani, Location on Net- 
works (M.LT. Press, Cambridge, MA, 1979). 

[9] S.L, Hakimi, Optimum locations of switching centers 
and the absolute centers and medians of a graph, 
Operations Res. 13 (1964) 462-475. 


116 O. Berman et al. | Network location 


{10] R.A. Howard, Dynamic Programming and Markov Pro- 
cesses (M.LT. Press, Cambridge, MA, 1960). 

{11] P. Jarvinen, I. Rajala and H. Sinervo, A branch-and- 
bound algorithm for seeking the p-median, Operations 
Res. 20 (1972) 173-178. 

[12] J.P. Jarvis, Optiniization in stochastic service systems 
with distinguishable servers, Technical Report IRP-TR- 
19-75, M.LT, Operations Research Center, Cambridge, 
MA (1975). 

{13] L. Kieinrock, Queuing systems, Vol. 1: Theory (Wiley, 
New York, 1975) Chapter 5S. 

{14] R.C. Larson, A hypercube-queuing model for facility 
location and redistricting in urban emergency services, 
Comput. Operations Res. 1 (1974) 67—95. 


[15] R.C. Larson, Structural system models for locational 
decisions: an example using the hypercube queuing 
model, in K.B. Haley, Ed., OR’78 (North-Holland, 
Arasterdam, 1979) 1054-1091. 

[16] P.B. Mirchandani, Analysis of stochastic networks in 
emergency service systems; Technical Report IRP-TR- 
15-75, M.LT. Operations Research Center, Cambridge, 
MA (1975). 

[17] P.B. Mirchandani and A.R. Odoni, Locations of medians 
on stochastic networks, Transportation Sci. 13 (1979) 
86-97. 

[18] C. ReVelle and R.W. Swain, Central facitities location, 
Geograph. Anal. 2 (1970) 30-42, 


